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ABSTRACT 

The  observation  of  the  motion  of  the  Earth's  spin  axis  with 
respect  to  the  crust  has  been  done  continuously  since  the  latter 
part  of  1899  by  the  International  Latitude  Service.  With  the  ad¬ 
vent  of  n.v  technologies,  new  determinations  of  polar  motion  have 
been  possible.  Doppler  tracking  of  the  Navy  Navigation  Satellites 
has  provided  estimates  of  the  polar  motion  on  a  permanent  basis 
since  1969.  Currently,  these  estimates  are  done  at  the  Defense 
Mapping  Agency  and  are  distributed  to  several  agenc'es  including 
the  Bureau  International  de  l'Heure  (BIH),  whi  •-  the  respon¬ 
sibility  of  centralizing  polar  motion  data.  The  ;  -CAR  Global 
Positioning  System  (GPS)  is  a  new  navigation  sate  111  .>  system  which 
will  eventually  replace  the  existing  Navy  Navigation  Satellite 
System.  As  a  byproduct  of  the  orbit  estimation  process  for  the 
GPS  satellites,  values  for  the  position  of  the  pole  are  determined. 
In  this  paper  tne  two  different  methods  for  computing  the  pole's 
position  from  satellite  observations  are  described.  The  most  re¬ 
cent  results  from  each  method  are  compared  to  each  other  and  to 
the  standard  BIH  values. 

\1msii2n 

Since  the  late  nineteenth  century,  it  has  been  recognized 
that  the  instantaneous  spin  axis  of  the  Earth  moves  with  respect 
to  the  geographic  pole  of  the  Earth's  crust.  This  "polar  motion" 
was  predicted  by  Euler  in  1752,  but  was  not  conclusively  observed 
until  Kiinstner's  work  in  1884-86.  The  International  Latitude  Ser¬ 
vice  (ILS)  was  established  in  1899  to  continuously  monitor  the  mo¬ 
tion  of  the  pole  by  making  systematic  determinations  of  latitude. 
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The  basic  idea  was  to  make  visual  zenith  telescope  observations  of 
stars  at  observatories  on  a  common  parallel  ( 39°8 1  N)  at  well-dis¬ 
tributed  longitudes  using  the  same  reduction  procedures.  The  ori¬ 
gin,  from  which  ILS  polar  motion  is  measured,  is  called  the  Con¬ 
ventional  International  Origin  (CIO).  The  CIO  pole,  which  is  de¬ 
fined  by  the  mean  latitude  of  the  five  ILS  observatories,  is  the 
mean  position  of  the  true  celestial  pole  from  1900  to  1905. 

With  the  development  of  new  astrometric  instruments  during 
the  1950s  and  their  use  at  the  national  observatories,  it  was  de¬ 
cided  that  latitude  measurements  made  at  stations  outside  the  ILS 
should  be  incorporated  for  polar  motion  determination.  Therefore, 
the  ILS  was  reorganized  into  the  International  Polar  Monitoring 
Service  (IPMS)  in  1962  with  the  charter  of  deriving  polar  motion 
from  latitude  and  universal  time  data  of  all  astronomical  instru¬ 
ments  . 

The  Bureau  International  de  l'Heure  (BIH) .was  created  in  1912 
to  unify  time  by  publishing  the  time  of  emission  (in  univeral  time) 
of  radio  time  signals.  With  the  advent  of  atomic  time  standards 
in  1955  and  the  need  to  account  for  polar  motion  in  a  timely  manner 
for  UT1  determination,  the  BIH  began  to  determine  its  own  set  of 
coordinates  of  the  pole  from  latitude  data  supplied  by  participa¬ 
ting  observatories.  Until  1972  the  BIH  used  only  the  same  data  as 
the  IPMS  so  a  complete  over Ian  of  function*  with  the  IPMS  existed. 
Since  1972,  pole  positions  baaed  oh  Doppler  observations  of  the 
Navy  Navigation  Satellites  have  been  incorporated  into  the  BIH 
global  solution.  A  detailed  description  of  the  method  used  by  the 
PIH  for  polar  motion  determination  was  given  by  Peissel  (1980). 

Research  at  the  Naval  Weapons  Laboratory  by  Anderle  and  Beu- 
glass  (1970)  demonstrated  that  it  was  possible  to  use  Doppler  ob¬ 
servation?  of  Navy  Navigation  Satellites  to  compute  pole  positions. 
Doppler  solutions  of  pole  position  have  been  distributed  by  the 
Dahlgren  Polar  Monitoring  Service  since  1969.  The  pole  positions 
are  a  byproduct  of  the  orbit  computation  process.  Hence,  when  the 
responsibility  of  computing  the  Navy  Navigation  Satellite  orbits 
was  transferred  to  the  Defense  Mapping  Agency  (DMA)  in  1975,  the 
derivation  of  pole  position  was  also  carried  out  by  DMA.  The  DMA*?-  - 
HTC  Polar  Monitoring  Service  (DPMS)  reports  are  distributed  by  tne 
Hydrographic/Topographic  Center  of  DMA  to  users  on  a  weekly  or 
monthly  basis  and  provided  to  the  U.  S.  Naval  Observatory  on  a 
weekly  basis  for  Inclusion  in  their  Time  Service  Announcement, 
Series  ?•  Additional  history  and  information  on  the  polar  motion 
services  has  been  given  by  Guinot  (1978). 

COMPUTATIONAL  METHODS  FOB  DOPPLER  DATA 

The  pole  determination  method  utilized  by  DMA  is  the  method 
adopted  by  the  Naval  Surface  Weapons  Center  (NSWC)  in  August  1971. 
The  brief  description  of  the  method  that  follows  is  taken  from  the 
detailed  description  of  the  observational  procedures , and  the  data 
reduction  techniques  given  by  Anderle  (1973). 


Doppler  observations  are  made  daily  by  a  network  of  approxi¬ 
mately  20  worldwide  tracking  stations  controlled  by  DMA  and  by  a 
network  of  4  U.  S.  tracking  stations  controlled  by  the  Navy  Astro¬ 
nautics  Group.  Sometimes  these  observations  are  supplemented  by 
tracking  data  from  portable  Doppler  receivers.  Two  types  of  Dop¬ 
pler  data  are  processed:  "sampled  Doppler  data"  which  gives  a  value 
of  the  Doppler  frequency  every  4  seconds  and  "continuous  count  in¬ 
tegrated  Doppler  (CC1D)  data"  which  gives  a  Value  of  range  differ¬ 
ence  each  10  to  20  seconds.  All  observations  taken  in  a  48-hour 
period  are  processed  by  the  CELEST  (O'Toole,  1976)  computer  pro¬ 
gram. 


Extensive  preprocessing  is  performed  or  the  Doppler  data. 

Data  are  rejected  if  transmission  errors  have  garbled,  sny  charac¬ 
ters  or  formats  or  if  unusual  time  gaps  appear  in  the  series  of 
observations.  Time  corrections  are  then  applied  to  the  data  and 
the  data  are  filtered  to  detect  gross  errors.  A  least-squares 
solution  for  a  modified  station  position  for  each  paws  is  performed 
to  linearize  residuals  for  filtering  purposes.  Residuals  for  the 
modified  station  position  are  rejected  if  they  are  excessive.  The 
process  continues  until  no  additional  points  are  rejected.  All 
passes  or  data  are  pnen  combined  and  passes  with  excessive  modified 
station  position  errors  are  rejected. 

After  the  preprocessing  of  two  day's  worth  of  data  is  done,  a 
least-squares  solution  is  obtained  which  includes  the  6  constants 
of  orbital  integration,  a  drag  scaling  factor  for  each  day,  a  fre¬ 
quency  and  a  tropospheric  refraction  scaling  factor  for  each  pass, 
the  2  components  of  the  pole  position,  and  the  coordinates  of  any 
mobile  observing  station.  The  least-squares  solution  is  based  on 
differences  between  the  observations  and  computed  data  which  cor¬ 
responds  to  a  predicted  satellite  orbit.  The  initial  conditions 
for  the  equations  of  motion  come  from  the  previous  day's  orbit  fit. 
The  integration  scheme  is  a  12th-order  Cowell  with  a  one-minute 
stepsize.  The  mathematical  model  includes  the  Earth's  gravitation¬ 
al  field,  atmospheric  drag,  solar  radiation  pressure,  luni-solar 
gravity  perturbations,  snd  solid  Earth  tidal  forces.  The  integra¬ 
tion  is  done  in  the  true-of-date  system. 

Two  satellites  are  currently  being  tracked  for  precise  ephem- 
eris  computations.  The  primary  satellite  is  1970-67A  and  the  sec¬ 
ondary  satellite  is  1967-92A.  Approximately  140  passes  of  data 
are  available  for  each  satellite  during  the  48-hour  time  span.  The 
current  accuracy  of  the  orbit  determination  as  measured  by  the  dis¬ 
continuity  in  the  ephemerls  at  the  ends  of  eaoh  48-hour  fit  and  the 
root  mean  square  of  the  residuals  of  the  least  squares  adjustment 
is  approximately  3  meters.  This  value  agrees  with  the  previous 
value  quoted  by  Bowman  and  Leroy  (1976). 

The  Initial  values  for  the  polar  motion  parameters  come  from  doing 
a  least  squares  fit  of  the  following  function  using  the  most  recent 
200  day's  worth  of  calculated  polar  motion  parameters. 


f(t)  =  A  +  E  cos(*>t)  +  C  sin(&>t) 
w  =  (21T/420)  where  420  =  Chandlerlan  period 


The  results  of  carrying  out  separate  fits  for  the  x-  and  y-compon- 
ents  of  polar  motion  are  used  to  predict  new  values  for  the  next 
seven  days . 


ACCURACY  OF  DOPPLER  POLE  POSITIONS 

Several  investigators  have  examined  the  error  sources  associ¬ 
ated  with  the  Doppler  results  and  the  accuracy  of  the  pole  positions. 
Anderle  (1973)  discussed  the  early  results  and  error  sources.  Oes- 
terwinter  (1979)  has  given  a  comprehensive  review  of  Doppler  results 
through  1977.  Bowman  and  Leroy  X1976)  hdve  performed  a  spectral 
analysis  of  the  x-  and  y-components  of  pole  position.  Guinot  (1979) 
has  examined  irregularities  of  the  polar  motion  by  comparing  clas¬ 
sical  astrometry  results  to  the  DMA  Doppler  results.  Gambis 
Nouel  (I960)  have  discussed  the  results  of  the  MEDOC  experiment 
which  was  an  independent  determination  of  pole  position  using 
Doppler  data.  Oesterwinter  gives  the  standard  deviation  of  the 
two-day  solution  with  respect  to  the  DPMS  five-day  mean  as  40  centi¬ 
meters,  and  the  standard' deviationof  the  five-day  mean  (standard- er¬ 
ror  ).  as  20  centimeters'.  Systematic  errors  are  due  to  an' inadequate 
knowledge  of  the  gravity  field  and  results  are  affected  by  atmos-*  • 
pfaeric  drag  variations  and  changes  in  the  station  network.  .  , 

The  accuracy  of  the  pole  positions  as  reflected  in  the  standard 
deviations  of  the  two-day  solutions  with  respect  to  the  five-day 
means  had  been  steadily  increasing  since  1967.  However,  the  results 
for  1981  show  a  much  larger  scatter.  The  average  values  for  the  x- 
and  y-components  were  84  (08027)  and  64  (0B021)  centimeters,  respec¬ 
tively.  The  larger  scatter  in  the  determinations  of  pole  position 
is  correlated  with  high  values  for  the  atmospheric  drag  parameters. 

NA7STAR  GLOBAL  POSITIONING  SYSTEM 

The  NAVSTAR  Global  Positioning  System  (GPS)  is  an  advanced  sat¬ 
ellite-based  navigation  system  currently  under  development  by  the 
Department  of  Defense.  It  will  ultimately  provide  position  and  ti¬ 
ming  information  accuracies  on  the  order  of  meters  and  tens  of  nan¬ 
oseconds  continuously  on  a  worldwide  basis  to  a  variety  of  users. 

The  GPS  consists  of  three  major  segments:  the  Space  Segment,  the 
Control  Segment,  and  the  User  Segment. 

The  Space  Segment  contains  the  satellites.  In  its  operational 
configuration  the  GPS  satellite  constellation  will  consist  of  18 
Navigation  development  Satellites  (NDSs)  in  circular,  12-hour  orbits 
(20,200  kilometers)  with  inclinations  of  55  degrees.  The  NDSs 
transmit  coded  signals  at  1575.42  and  1227.6  MHz  which  consist  of  a 
precise  navigation  signal,  a  course  navigation  signal,  and  data  such 
as  satellite  ephemerldes,  satellite  clock  bias  information,  andat- 
mospherlc  propagation  correction  data.  Currently,  the  useful  GPS 


satellite  constellation  consists  of 
hour  orbits  with  inclinations  of  63 

The  Control  Segment  contains  four  monitor  stations,  an  upload 
station,  and  a  Master  Control  Station  (MCS).  It  is  responsible  for 
tracking  the  satellites,  preparing  the  navigation  data  for  the  sat¬ 
ellites,  and  uploading  the  data  into  a  satellite  navigation  pro¬ 
cessor.  This  process  Is  done  at  least  once  per  day. 

The  User  Segment  contains  the  various  user  equipment  sets  which 
demodulate  the  navigation  data  sent  by  the  satellites  and  process 
them  to  determine  the  user's  three-dimensional  position,  velocity , 
and  system  time.  Examples  of  typical  user  performance  are. given  by 
Henderson  and  Strada  (1980).  Additional  information  on  the  princi¬ 
ple  of  operation  and  system  characteristics  of  GPS  are  given  by 
Milliken  and  Zoller  (1978). 

The  GPS  program  is  currently  in  the  full-scale  development  and 
system  test  phase,  in  two  years  the  decision  whether  to  enter  the 
production  and  development  phase  leading  to  full  GPS  capability  will 
be  made.  Because  the  program  is  in  the  test  phase,  it  seems  appro¬ 
priate  to  examine  whether  useful  information  can  be  obtained  for 
pole  position  as  a  byproduct  of  the  orbit  computation  process. 


The  MCS  at  Vandenberg  APB,  California,  completely  controls  the 
operation  of  the  Control  Segment  of  the  GPS.  The  primary  function 
of  the  MCS  is  to  provide  the  NDSs  with  precise  navigation  data. 
Tracking  data  consisting  of  pseudoranges  (true  slant  ranges  from 
the  user  to  the  satellite  plus  propagation  delays  and  time  biases) 
and  delta  pseudoranges  (integrated  Doppler)  are  collected  by  the 
monitor  stations  and  sent  to  the  MCS  for  processing.  The  data  are 
preprocessed,  smoothed,  and  then  run  through  a  Kalman  filter  to  es¬ 
timate  the  NDS  ephemerides  and  clock  behavior  so  that  predicted- eph- 
emerides  can  be  uploaded  to  the  satellites. 

The  MCS  ephemeris  determination  processing  is  summarized  func¬ 
tionally  in  Fig.  1.  The  monitor  stations  collect  the  following  data 
for  each  satellite  in  view  (nominally  every  6  seconds):  pseudo¬ 

range  measurement,  L1-L2  pseudorange  difference  measurement,  and 
delta  pseudorange  (integrated  Doppler)  on  Li.  In  addition  to  the 
measurement  data,  meteorological  and  satellite  (SV)  navigation  data 
are  sent  to  the  MCS.  The  preprocessor  (FRE?)  unpacks  and  scales  the 
measurement  data,  edits  bad  data,  corrects  the  measurement  time  tag, 
and  corrects  the  measurements  for  antenna  phase  center  offsets,  ion¬ 
ospheric  delays,  tronosnheric  effects,  relativistic  effects,  and 
Earth  rotation  during  signal  transit.  The  smoother  (SMOOTH)  proc- . 
esses  all  measurements  contained  within  the  15  minute  Kalman  filter 
cycle  for  each  satellite-monitor  station  pair.  Wild  points  are  ed- 
itted.  The  remaining  pseudorange  and  delta  pseudorange  measure¬ 
ments  are  fit  by  least  squares  to  a  polynomial  and  1  smoothed  pseudo¬ 
range  residual  and  1  smoothed  delta  pseudorange  residual  with 


I  satellites  in  circular,  12- 
degrees  . 
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appropriate  statistics  are  sent  on  to  the  Kalman  filter  (KALMAN). 

The  Kalman  filter  processes  the  smoothed  measurements  to  estimate 
NDS  position  and  velocity  residuals,  NDS  solar  radiation  pressure 
model  scaling  parameter  residuals,  NDS  clock  offset  states,  monitor, 
station  clock  offset  states,  tropospheric  delay  residual  states, 
and  polar  wander  residual  states.  KALMAN  is  a  square-root,  sequen¬ 
tial  Kalman  filter  process.  The  predictor  (PHEDICT)  uses  the  most 
"ecent  Kalman  estimates  of  the  residuals  to  predict  the  Earth-fixed 
Cartesian  coordinates  of  the  satellite  position  which  can  then  be 
uploaded  to  the  satellite.  The  Kalman  filter  implementation  re¬ 
quires  a  reference  trajectory  which  can  be  generated  by  an  integra¬ 
tor  (INTEGRATE)  or  by  the  Naval  Surface  Weapons  Center  (NSWC).  The 
relevant  portions  of  the  reference  trajectory  needed  for  PHEP,  KAL¬ 
MAN,  and  PHEDICT  are  obtained  by  the  extractor  (EXTRACT) . 

The  reference  trajectory  provides  the  initial  estimate  of  the 
satellite  trajectory  about  which  the  perturbations  are  computed  by 
the  Kalman  filter  and  the  partial  derivatives  which  are  also  needed 
by  the  Kalman  filter.  The  smoothed  tracking  data  are  transmitted 
to  NSWC  each  day.  A  least-squares  batch  fit  using  range  difference 
data  is  used  to  generate  the  reference  ephemerides  every  two  weeks 
at  NSWC.  Predicted  values  of  the  pole  position  are  also  included 
on  the  reference  tape.  These  values  are  obtained  by  fitting  a 
function  similar  to  the  Doppler  function  but  also  including  a  sina 
and  a  cosine  term  with  an  annual  period.  A  Chandler  period  of  435 
days  is  used  and  the  most  recent  435  day's  worth  of  BIH  data  are 
used  for  the  fit  (Carr,  1981).  For  the  NDSs  on  atomic  standards, 
typical  errors  for  a  two-week  prediction  versus  the  fit  span  are  5 
meters  root-mean-square  (rms)  in  the  radial  and  normal  directions 
and  50  meters  rms  in  the  tangential  direction. for  the  reference 
trajectories. 

When  the  Kalman  filter  was  designed  at  the  MCS,  it  was  decided 
to  not  describe  the  dynamics  of  the  polar  wander  states.  Therefore, 
the  states  are  constant  and  propagate  in  time  via  the  identity  ma¬ 
trix.  The  Kalman  estimates  oi  tnese  states  (based  on  ranging  mea¬ 
surements)  wiijL  unange  wltn  time,  but  epoch  updates  (updating  the 
reference  values)  will  not  change  the  total  state.  The  following 
equation  describes  the  epoch  update  process  for  polar  wander. 

PWCn  *  PW0  +  PWco  -  Wn 

where  FWo  =  old  reference  (NSWC)  value  of  polar  wander 
PWn  3  new  reference  (NSWC)  value  of  polar  wander 

PWco  3  current  estimate  of  polar  wander  residual  state 
relative  to  the  old  reference  value 

PWcn  3  current  estimate  of  polar  wander  residual  state 
relative  to  the  new  reference  value 

The  first  two  terms  on  the  right-hand  side  of  the  equation  represent 


the  total  polar  wander  state.  The  process  noise  matrix  describes 
the  estimated  uncertainty  and  correlations  of  the  Kalman  filter 
states  due  to  unmodeled  effects  or  modeling  errors.  The  process 
noise  matrix  for  the  polar  wander  states  assumes  that  the  residual 
states  behave  like  stationary  independent  white  noise  random  pro¬ 
cesses  over  the  periods  between  successive  NSWC  reference  values 
of  polar  wander.  Nominal  values  for  the  process  noise  correspond 
to  a  position  uncertainty  of  0.?  meter  over  a  period  of  a  day. 

The  values  of  the  pole  position  based  on  MCS  Kalman  filter 
processing  were  obtained  from  a  dump  of  the  Kalman  estimate  file 
just  before  the  first  scheduled  uploads  of  the  NDSs  each  day.  As 
long  as  the  current  estimate  of  the  polar  wander  residual  state  is 
not  reset  to  zero,  the  total  polar  wander  state  is  independent  of 
the  reference  value  supplied  by  NSWC. 

Additional  information  on  the  operation  of  the  MCS  and  the 
Control  Segment  has  been  given  by  Russell  and  Schaibly  (1978)  and 
Varnum  (1982). 


ACCURACY  OP  GPS  POLS  POSITIONS 

Anderle ,  Beuglass,  and  Carr  (1981)  have  computed  pole  position 
and  Earth  rotation  for  six  7-day  spans  of  CPS  tracking  data  observed 
during  the  MERIT  short  campaign.  They  concluded  that  both  standard 
errors  and  consistency  of  the  results  were  somewhat  worse  than 
those  obtained  from  the  Navy  Navigation  Satellite  data.  They  give 
standard  deviations  of  the  x-  and  y-components  ranging  from  26  to 
85  centimeters  depending  on  the  assumed  uncertainty  of  station  bias 
and  frequency. 

The  standard  deviations  of  the  daily  solutions  with  respect  to 
the  five-day  means  are  quite  small.  The  average  value  for  both  the 
x-  and  the  y-comnonent  is  6  (0!002)  centimeters.  These  small  values 
are  indicative  of  the  constraints  on  the  solution  in  the  MCS  Kalman 
filter. 


COMPARISON  OP  RESULTS 

The  pole  paths  computed  by  Doppler  observations  of  the  Navy 
Navigation  Satellites  (DMA),  by  MCS  Kalman  filter  processing  of  NDS 
data  (MCS),  and  by  the  BIH  are  shown  in  Pig.  2.  The  DMA  data, , the 
MCS  data,  and  the  BIH  data  are  plotted  as  dots,  pluses,  and  circles, 
respectively.  The  scatter  in  the  DMA  data  with  respect  to  the  BIH 
is  quite  large.  In  addition  the  DMA  data  appear  to  be  systematic¬ 
ally  deviating  from  the  BIH  in  the  final  part  of  the  year.  The  MCS 
data  show  a  signif'cant  biai  with  respect  to  the  BIH  and  the  DMA 
data  for  most  of  t*  -  vear  The  BIH  data  are  the  Circular  D  final 
values.  Graphs  of  on*  cc._ponents  of  pole  position  as  a  function  of 
time  are  shown  in  Figures  3  and  The  MCS  data  fit  the  BIH  values 
better  than  DMA  data  for  the  x-coraponent  while  just  the  opposite  is 
true  for  the  y-component.  Residuals  of  the  different  estimates  of 
pole  position  are  given  in  Figures  5  and  6.  Table  1  presents  a 
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summary  of  the  root-mean-square  differences  of  the  residuals  given 
in  Figures  5  and  6. 


Pole  Position  Differences  (arc  sec) 


MCS 

-  BIH 

DMA 

-  BIH 

MCS 

-  DMA 

X 

Y 

JL_ 

_ _ Y 

X 

Y 

mean 

difference 

.009 

•  035 

-.018 

i 

e 

o 

o 

.027 

• 

o 

00 

std. 

deviation 

.005 

.012 

.015 

.012 

.016 

.013 

std. 

error 

.001 

.002 

• 

o 

o 

.002 

.003 

.002 

TABLE  1.  COMPARISON  OF  MEAN  ESTIMATES  OF  COORDINATES  OF  THE  POLE 

SfflHBABX 

Pole  positions  have  been  computed  as  a  byproduct  of  the  orbit 
estimation  process  for  the  GPS  satellites  at  the  Master  Control 
Station  (MCS).  The  precision  of  the  daily  values  indicates  an  un¬ 
certainty  of  05002  or  6  centimeters.  However,  systematic  errors 
are  apparent  in  the  results.  A  bias  of  about  1  meter  exists  between 
the  MCS  pole  positions  and  the  positions  given  by  the  Bureau  Inter¬ 
national  de  l'Heure  (BIH)  for  the  y-component.  When  the  MCS  re¬ 
sults  are  compared  to  the  BIH  values,  the  scatter  about  the  mean 
difference  is  much  smaller  than  for  the  DMA  pole  positions.  In  ad¬ 
dition  to  the  systematic  errors,  some  periodicities  mav  exist  in 
the  data  but  the  data  span  of  MCS  values  is  still  too  short  to  do 
a  comprehensive  analysis. 

Preliminary  analysis  of  MCS  pole  positions  has  given  encourag¬ 
ing  results.  However,  more  data  need  to  be  analyzed  and  additional 
studies  need  to  be  completed  before  definitive  statements  can  be 
made  about  the  causes  of  the  biases  and  the  appropriateness  of  in¬ 
cluding  these  data  in  a  global  solution  for  polar  motion. 
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